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Abstract. We investigate the impact that warm dark matter (WDM) has in terms of 21cm intensity 
mapping in the post-reionization Universe at z = 3 — 5. We perform hydrodynamic simulations for 
5 different models: cold dark matter and WDM with 1,2,3,4 keV (thermal relic) mass and assign the 
neutral hydrogen a-posteriori using two different methods that both reproduce observations in terms 
of column density distribution function of neutral hydrogen systems. Contrary to naive expectations, 
the suppression of power present in the linear and non-linear matter power spectra, results in an 
increase of power in terms of neutral hydrogen and 21cm power spectra. This is due to the fact that 
there is a lack of small mass halos in WDM models with respect to cold dark matter: in order to 
distribute a total amount of neutral hydrogen within the two cosmological models, a larger quantity 
has to be placed in the most massive halos, that are more biased compared to the cold dark matter 
cosmology. We quantify this effect and address significance for the telescope SKAl-LOW, including 
a realistic noise modeling. The results indicate that we will be able to rule out a 4 keV WDM model 
with 5000 hours of observations at z > 3, with a statistical significance of > 3 cr, while a smaller 
mass of 3 keV, comparable fo presenf day consfrainfs, can be ruled ouf af more fhan 2 a confidence 
level wifh 1000 hours of observations af z > 5. 
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1 Introduction 

The ACDM cosmological model has been shown to be very successful in explaining a wide variety 
of observables: the anisotropies in the cosmic microwave background, the clustering of galaxies at 
low redshift, the abundance of galaxy clusters, the baryonic acoustic oscillations observed both in the 
spatial distribution of galaxies and in the Lyman-a forest. In this model, the dark matter is assumed 
to be cold (CDM), i.e. with negligible thermal velocities on all scales at high redshift. Although some 
tensions arise when combining different scales, a general picture emerge for which the simplest six 
parameter model provides a good fit to all these data sets [1]. 

On small scales however, some predictions of the ACDM model seems to be at odds with 
observations. Among others, two major problems show up: the so-called core-cusp problem and the 
subhalos abundance problem. The former arises because N-body simulations predicts the density 
profile of dark matter halos to be cuspy [2, 3], whereas observations find them, in some galaxies, 
cored [4-8]. The second problem comes because N-body simulations predicts, for Milky Way halos, 
the existence of a larger number of subhalos than the number of galaxy satellites found in the vicinity 
of our galaxy. More recently, this apparent crisis of the ACDM scenario has been casted in terms 
of the ’’too-big-to-fail” problem: namely the dynamical properties of the most massive subhalos of 
a Milky Way simulated halo are at odds with those of the observed satellites (dwarf galaxies), there 
must be a mechanism that prevent the star formation in this massive halos [9]. The most natural 
explanation would probably be to invoke some yet to be understood baryonic physical processes that 
deeply affect the properties at small scales [10-12]. 

However, there is also a possible way to alleviate at least some small scales problems, while 
preserving the success of CDM on large scales, by assuming that the dark matter particles have 
intrinsic thermal velocities (although there are works that disagree, see for example [13, 14]). On 
one hand, the presence of thermal velocities will avoid the cuspy clustering of dark matter particles 
within gravitational potential wells, allowing the formation of cores [15-19] in the density profile. 
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although most of the works have shown that an unrealistic low value of the WDM mass is needed 
to match observations. On the other hand, the collapse of dark matter halos can only take place 
for perturbations larger than the typical free-streaming length of the dark matter particles, thus, the 
presence of thermal velocities in the dark matter particles will naturally suppress the abundance of 
low mass subhalos [20-22] . 

While structure formation in a universe dominated by dark matter particles with large thermal 
velocities (hot dark matter) will proceed from top to bottom, the bottom-top paradigm is preserved 
in the case of dark matter particles with relatively low thermal velocities that could be considered 
as warm dark matter (WDM). WDM candidates comprise sterile neutrinos and gravitinos or other 
thermal relics. 

Given the fact that WDM can not substantially cluster on scales smaller than its free streaming 
scale, the matter power spectrum in a cosmological model with WDM will exhibit, on small scales, 
a cut-off. Because the Lyman-a forest can be used to probe the shape and amplitude of the matter 
power spectrum on small scales, the current tightest constrains on the WDM mass come from such 
observations: m-yvoM ^ 3.3 keV(2(T) [23]. 

Nowadays, our constraints on the matter power spectrum arise from different probes such as 
the galaxy clustering, weak leasing, abundance of galaxy clusters and the Lyman-a forest. In the 
near future, observations of the redshifted 21cm radiation from neutral hydrogen will allow us to 
constrain the matter power spectrum on redshifts not accessible with the above observables with an 
unprecedented precision [24, 25]. The idea is to study the statistical properties in the integrated 
emission from unresolved galaxies using low angular resolution surveys; this technique is called 
intensity mapping. This observable is promising: it is at high redshift and thereby should be closer 
to sample the linear cut-off in the matter power spectrum induced by WDM; the scales and redshifts 
addressed here should only be marginally affected by baryonic processes like feedback; the total 
amount of neutral hydrogen is constrained by data, i.e. absorption lines. 

In order to extract the maximum information from intensity mapping surveys, it is critical to 
model the spatial distribution of neutral hydrogen. In this paper we investigate the differences in the 
spatial distribution of neutral hydrogen between models with cold and warm dark matter. Since the 
quantity directly probed by radio-telescopes is the 21cm power spectrum, we study the signatures 
imprinted by WDM on this observable. We do so by running high resolution hydrodynamic simu¬ 
lations. We also compute the sensitivity with which the future Square Kilometre array ^ (SKA) will 
constrain the 21cm power spectrum, and therefore, how much it will be able to discriminate among 
different cosmological models. 

The scientific motivation behind this work is actually more general than the WDM framework 
we are investigating since we are aiming at providing a quantitative interpretation of the clustering 
signal in intensity mapping in a completely new regime of scales and redshifts, with very little (if 
any) overlap with other observational probes. 

We notice that the signatures left by WDM on the 21cm power spectrum has already been stud¬ 
ied in [26], which focused their study in the reionization era, with its richer astrophysical processes 
while in this work we investigate the post-reionization era. Furthermore, the impact of WDM on the 
higher redshift regime has been recently explored by [27, 28], and more exotic annihilating DM and 
decaying WDM are taken into consideration by [29] . 

This paper is organized as follows. In Sec. 2 we present the hydrodynamic simulations used 
in the present work. In this section we also study the differences in the spatial distribution of matter 
between the models with CDM and WDM. The two different methods used to simulate the spatial 

* https ://www. skatelescope.org/ 
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distribution of neutral hydrogen are described in Sec. 3, while the spatial distribution of neutral 
hydrogen in the models with CDM and WDM is investigated in Sec. 4. The capability of the future 
SKAl-LOW array to distinguish between the different models is fully quantified in Sec. 4. Finally, a 
summary of the work carried out in this paper and the main conclusions of it are presented in Sec. 5. 

2 Simulations and matter distribution 

In this section, we present first the hydrodynamic simulations in subsection 2.1, whereas the impact 
of WDM on the spatial distribution of matter, studied in terms of the matter power spectrum and the 
halo mass function is discussed in subsection 2.2. 

2.1 Hydrodynamic simulations 

Our simulation suite comprises a set of 5 high-resolution hydrodynamic N-body simulations run using 
the TreePM-tSPH code GADGET-III [30]. We follow the evolution of 512^ CDMAVDM and 512^ 
baryon (gas-i-stars) particles within simulation boxes of comoving sizes equal to 30 /i“^Mpc. The 
mass resolution of our simulations is ttidm = 1-50 x 10^ and mb = 2.74 x 10®The 

gravitational softening is set to 1/40 of the mean inter-particle linear spacing, i.e. 1.5 h~^kpc. We 
have simulated five differenf cosmological models: 1 model wifh CDM and 4 models wifh WDM, 
each of fhem wifh a differenf Fermi-Dirac mass of fhe WDM particles: 1 keV, 2 keV, 3 keV and 
4 keV (fhermal relic mass, see below). The values of fhe cosmological parameters are fhe same in 
all fhe models and in agreemenf wifh fhe lafesf resulfs of fhe Planck satellite [1]: Dm = 0.3175, 
Db = 0.049, Da = 0.6825, h = 0.6711, = 0.9624 and cjg = 0.834. We nofice fhaf fhe resolufion 

on our simulations is high enough fo provide converged resulfs when fhe HI is assigned eifher wifh 
fhe halo or particle based mefhod (see [31]). 

Sfar formation is modeled using fhe effecfive mulfi-phase model of Springel & Hemquisf [32] . 
The code also simulafes radiative cooling by hydrogen and helium and heating by a uniform Ulfra 
Violef (UV) background. Bofh fhe cooling roufine and fhe UV background have been modified fo 
obfain a desired fhermal history, corresponding fo fhe reference model of [33], which is kepi fixed in 
fhe following. 

The initial condifions were generafed af 2 : = 99 using fhe Zel’dovich approximation. For fhe 
model wifh CDM fhe mailer power speclrum and fhe Iransfer funcfions were computed using GAME 
[34] whereas for fhe models wifh WDM fhe power speclra, Pa{k), were calculaled as [35] 

( 2 . 1 ) 

where a sfands for eifher baryonic and fhe non-baryonic mailer componenls and wifh 

Tim(fc) = (l + (aA:)2-)-^/" (2.2) 

where v = 1.12 and 

= /.-■Mpo . (2.3) 

V"1WDM/ V 0-25 J J 

Furfhermore, for fhe cosmological models wifh WDM we have added fhermal velocifies, on lop of 
fhe peculiar velocifies, to fhe non-baryonic particles. For a given WDM particle, fhe modulus of fhe 
fhermal velocily veclor is drawn randomly from a Fermi-Dirac dislribulion wifh a mean equal fo [20] 

Vwdm( 4) = 0.012(1 + b ^ ' km/B (2.4) 
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Figure 1. Dimensionless power spectrum of total matter (DM + baryons), A^(fc) = P^{k) /{2Tr^), at 
redshifts z = 3 (left), z = A (middle) and z = 5 (right). In the top panels we plot with a solid black line the 
matter power spectrum of the model with CDM and with colored lines the results for the WDM models: 1 keV 
(dash-dotted blue), 2 keV (dashed red), 3 keV (dotted green) and 4 keV (dash-dotted magenta). The bottom 
panels show the relative differences between the WDM and the CDM models. We set the k range up to the 
Nyquist frequency (fc ~ 107 /iMpc”^). 

and the direction of the velocity vector is taken randomly. We analyzed snapshots at z = 3,3.5,4,4.5, 5 
and identified the dark matter halos and subhalos using the Friends-of-Friends (FoF) [36] and SUB- 
FIND [37, 38] algorithms, respectively. 

2.2 Impact of WDM on the matter distribution 

We now examine the impact of WDM on the spatial distribution of matter. In particular, we focus our 
attention on the matter power spectrum and the halo mass function. This high redshift regime z = 

3 — 5 is relatively new for probing WDM models in the sense that it partially overlaps with Lyman- 
a forest data, while is lower compared to the reionization and early structure formation regimes 
investigated by [26-28]. 

In the top panels of figure 1 we show fhe total matter power spectrum, P^{k), for each of our 
five different cosmologies. The bottom panels represent instead the relative differences between the 
WDM and the CDM models. Since WDM can not substantially cluster on scales smaller than its 
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Figure 2. Halo mass function (using the FoF catalogue) at redshifts z = 3 (left), z = 4 (middle) and z = 5 
(right). In the top panels the halo mass function for the CDM model is shown with a solid black line whereas the 
colored lines represent the results for the WDM models: 1 keV mass (dash-dotted blue), 2 keV (dashed red), 3 
keV (dotted green) and 4 keV (dash-dotted magenta). The bottom panels show the relative differences between 
the WDM and the CDM models. Error bars represent the uncertainty in the halo mass function assuming that 
the number of halos follows a Poissonian distribution. 


free-streaming length, the matter power spectrum in WDM cosmologies presents a relatively sharp 
cut-off on small scales; as expected, this suppression, while being less pronounced that the linear 
cut-off, is larger at higher redshift and for smaller WDM masses. The clustering on scales larger than 
the free-streaming length is not affected by the thermal velocities of the WDM particles, and thus, 
the amplitude of the matter power spectrum is the same for all the models. 

In figure 2 we show the halo mass function for the different models investigated in this paper. 
In the top panels we plot the halo mass function at 2 ; = 3 (left), z = 4 (middle) and z = 5 (right). The 
error bars represent the errors on the halo mass function assuming that the number density of halos 
follows a Poissonian distribution. In the bottom panel, we plot the relative difference between the 
WDM models and the model with CDM. As it can be seen in the figure, WDM induces a suppression 
on fhe abundance of low-mass dark mailer halos, wilh respecf fo Ihe abundance of halos in Ihe CDM 
model, which increases wilh redshift and decreases wilh Ihe mass of Ihe WDM parlicles. As we 
shall see laler, Ihis effecl will drive Ihe changes we observe in Ihe 21cm power speclrum belween Ihe 
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models with CDM and WDM. 

In deriving the mass function of WDM simulations, spurious fragmentation could produce an 
increase of halos at a low mass scale, see for example [22]. Our simulations, except the model with 
1 keV, are not affected by such effect, both because we artificially boost late time thermal velocities 
at the beginning of our simulations [39] and because we do not reach enough mass resolution for our 
WDM models, which have higher masses compared to those in the literature. This can be directly 
seen from the behaviour of the curves in figure 2, fhaf do nof show any characferisfics sfeep power law 
upfum for low halo masses. Moreover, as shown in [39], fhe mass function of fhe artificial clumps 
sfrongly depends on fhe resolution of fhe simulations; we checked fhe mass funclions resulfing from 
simulations of 1/8 fhe resolution we have (i.e. comoving box size of 60 /i“^Mpc) and we have 
verified fhe convergence for small halo masses. 

Overall, for WDM masses befween 3 and 4 keV we observe a suppression in power up fo 
~ 10% on very small scales and a reducfion in fhe number of 10® h~^MQ halos of fhe order of 
20-40% compared fo fhe CDM case. 

We til fhe WDM mass funclions obfained from fhe simulafions (FoF calalogue) by paramelriz- 
ing fheir deviafion from fhe mass function predicled by fhe Shefh & Tormen [40] formula^ 


rl/WDMf/Vf') = ^^WDM 
d log M ^ d log M 


yivi) 


(2.5) 


where a and Mq are free paramefers. The besf til values for fhe differenl WDM models and redshifls 
are shown in fhe fable 1. We find lhal fhe mass funclion of fhe 1 keV WDM cosmology is fhe only one 
fhaf is under predicled by fhe Shefh & Tormen model. We exclude fhe possibilily of fragmenfalion, 
because lacks fhe associafed sleeper raise for low halo masses (as already discussed above), 

allhough Ihis effecl could poinf fo a Iransifional and physical regime before fragmenfalion lakes place 
(see figure B1 in [39]). We also find lhal fhe mass function of fhe CDM model is over predicled by 
fhe Shefh & Tormen formula, alfhough only in fhe low mass end; fhe besf til values obfained using 
equalion 2.5 are shown in fhe fable 1 wilh cx) as fhe value of fhe WDM mass. 


3 Modeling the neutral hydrogen distrihution 

In Ihis secfion we describe fhe Iwo differenl approaches we use fo model fhe spalial dislribulion of 
neufral hydrogen. In subsecfion 3.1 we presenl fhe halo based melhod, which relies on fwo assump- 
fions: no HI is found oufside dark mailer halos and fhe amounl of neufral hydrogen inside fhe halos 
is a funclion of fheir lolal mass only. The so-called particle based melhod is summarized in subsec¬ 
fion 3.2: according fo Ihis modeling, HI is assigned fo every single gas particle in fhe simulafion, 
depending on fhe physical properties of fhe parficles Ihemselves. Bolh models have been extensively 
invesfigafed for inlensily mapping al similar redshifls in [31] and are meanl fo brackel conservalively 
fhe possible ways of populating halos wilh HI. 

3.1 Halo based method 

Here we follow one of fhe models described in Bagla ef al. [41]. Slarling from fhe snapshols of our 
N-body simulations, we first single out the dark matter halos and all the particles belonging to them, 
as identified by the Friends-of-Friends algorithm (FoF catalogue). We then assume that the amount 
of HI hosted by a dark matter halo of mass M is given by the function Mhi(M), i.e. we neglect any 

^This is achieved by computing the Sheth & Tormen mass function using the linear power spectrum of the WDM model. 
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tTTWDM [keV] 

z 

a 

Mo [10^ 

1 

3 

-1.72 

0.887 

4 

-1.92 

0.960 

5 

-2.06 

0.737 

2 

3 

0.571 

3.98 

4 

0.389 

3.80 

5 

0.286 

4.81 

3 

3 

0.830 

1.86 

4 

0.726 

1.28 

5 

0.427 

1.74 

4 

3 

0.958 

1.44 

4 

1.04 

0.890 

5 

0.522 

1.05 

oo 

3 

1.26 

1.04 

4 

2.77 

0.491 

5 

1.69 

0.455 


Table 1. Parameter values for the mass function fit of equation 2.5. 



Figure 3. The function Mhi(M) relating the HI mass and total mass of a dark matter halo of mass M, 
employed in the halo based model as described by the equation 3.1, for the CDM (dashed lines) and 1 keV 
WDM (dotted lines) cases, for redshifts z = 3,4, 5 (blue, red and green). The values of the parameter /a are 
0.0183, 0.0239, 0.0337 for the CDM, redshift 3 to 5, and 0.0325, 0.0491, 0.0835 for the 1 keV, redshift 3 to 5. 


environmental effects. We then split the total HI mass of the dark matter halo equally among all its 
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gas particles. Following the model 3 of Bagla et ah, we model the function Mhi(M) as: 


Mhi(M) 


f M 

•13 1 + (M/Mn,ax) 

0 


if Mmin < M 
Otherwise. 


(3.1) 


The values of the free parameters M^ax and Mm^n are chosen in correspondence dark matter halo 
circular velocities Vcuc of 200 and 30 km respectively, using the relation 


M = 10^°Mo 


Vc\r 


60 kms 


-1 


n + z 

V 4 


- 1.5 


(3.2) 


derived from virialisation arguments. The lower mass cut-off Mmin takes into account that a mini¬ 
mum hydrogen density (clustered in a minumum halo potential well) is needed to have self-shielding 
and prevent the gas to be fully ionised. The fall-off above M^ax is introduced to mimic observations: 
galaxies in dense environments as galaxy clusters are Hl-poor. For further details on this model we 
refer the reader to [41]. The magnitude of Mm\n is of the order of 10® so to distribute 

HI consistently, our simulations need to resolve dark matter halos up to this mass halos or lower. 
We achieve this by setting the box size of our simulations to 30 /i“^Mpc, among the highest values 
possible considering the number of particles, 2 x 5123, we are dealing with. 

We are left with a single free parameter, /a, whose value is chosen by demanding that the total 
amount of HI in our boxes reproduces the observational constraints arising from the abundance of 
DLAs at high redshift: Hhi(' 2) = 10“3, i.e. we impose: 


M- 

^3 E 1-, (3-3) 

-^max 

where L = 30 h^^Mpc is the simulation box size, p® = 2.78 x IO^^/i^MqMpc"® is the present 
day critical density of the universe, 0(x) is the Heaviside step function and the index i run over all 
the dark matter halos of the simulation. We stress that by defining /s in this way for each snapshot 
and simulation, we assure to have the same reference value of Hhi independently of the redshift and 
cosmology we analyze, within the halo based method. We will show that this is not the case in the 
particle based method. 

We show the function Mhi(M) that the halo based method prescribes for the models with 
CDM and 1 keV WDM in figure 3. The Mhi(M) funcfion sfeeply increases up fo Mmax> then 
displays a plafeau. The values of and Mm ax do nof change among differenf cosmologies, buf 
fhe normalizafion (fhrough fhe paramefer /a) does: if gefs higher for warmer DM scenarios. 


3.2 Particle based method 

By employing this method we assign neutral hydrogen to all gas particles in the simulation according 
to their physical properties. Compared to the halo based method, we do not rely on any definition of 
dark matter halo or on any assumption on the amount of HI outside halos and therefore we can also 
predict the amount of HI in places different to halos as filaments and cosmic voids. 

We proceed as depicted in [42]: for every single gas particle in the simulation we compute 
the neutral hydrogen fraction in photo-ionization equilibrium with the Ultra Violet (UV) background 
and correct that fraction to account for both HI self-shielding and formation of molecular hydrogen. 
Hydrogen of star forming particles is split in totally neutral for the cold phase, and totally ionized for 
the hot phase, i.e. we consider than the HFH fraction of star forming particles in photo-ionization 
equilibrium is equal to its multi-phase cold gas fraction. The strength of the UV background, at a 







particular redshift, is corrected by demanding that the mean flux of the Lyman-a forest reproduces 
the observations [43], with the approach described in [31]. For every gas particle we then compute 
the radial column density profile Niii{r) making use of the SPH spline kernel W(R): 


Nm{r) 



(3.4) 


with m/mH being the gas particle mass in units of hydrogen atom mass, HI/H is the neutral hydrogen 
fraction in photo-ionization equilibrium and rspH the SPH smoothing length. If exists a radius r such 
that Niii{r) = cm“^, then the sphere from i? = 0 to ii = r is considered to be self-shielded 

agains the external UV radiation and its HFH fraction is set to 0.9. For the surrounding spherical 
shell {R > r) it still holds the neutral fraction HI/H coming from the photo-ionization equilibrium 
assumption. 

The last step consists in correcting the HI/H fraction computed in the previous step to account 
for the formation of molecular hydrogen. Using the density and internal energy, we compute the 
pressure of gas particles P and use it to correct the HI/H fraction computed previously using the 
observed relation between the surface densities of HI and H 2 with the pressure of the disk galaxy 
[44]: 



with k-Q the Boltzmann constant. H 2 is assigned only to star forming particles. We note that we are 
implicitly assuming that the above relationship, which arises from observations at low-redshift, holds 
also at high redshift which may not be the case (see for instance [42, 45]). 

4 Results 

In this section we use the methods described above to model the spatial distribution of neutral hydro¬ 
gen in our five different cosmological models. 

An overall picture of the spatial distribution of matter and HI analyzed in this work is given in 
figure 4, where we confront the CDM scenario (left column) with the 1 keV WDM (right column) at 
redshift z = 3. The difference among the two cosmologies can be noted in the total matter distribution 
(top panels), where the CDM displays much more clustered structures while the spatial distribution 
of HI also differ in the two cosmologies (middle and bottom panels). By using the halo based, HI 
is present just in halos (middle panels), while by employing the particle based (bottom panels), the 
HI distribution is smoother since HI is assigned to every single gas particle in the simulation. The 
differences in the neutral hydrogen clustering properties can also be noticed by visual inspection. 

We perform a more quantitative analysis of the HI distribution. First, in the subsection 4.1, we 
compare the HI column density distribution function that we obtain, for each cosmological model, 
with observations. In the subsections 4.2 and 4.3 we investigate the differences in terms of HI power 
spectrum and of the HI bias ~ PHi{k)/Pm{k)- The 21cm power spectrum is computed in the 

subsection 4.4, where we also investigate the accuracy with which the SKA 1-LOW telescope will be 
able to measure it. 

4.1 HI column density distribution function 

Up to present days our understanding of the morphology of the hydrogen content in the high redshift 
universe comes mainly from quasar absorption spectra, where the presence of the Lyman-a transition 
lines is a conclusive evidence of intervening gas between us and the sources. From these observations 
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Figure 4. Spatial distribution, at redshift 2 = 3, of the density contrast of total matter (DM + baryons) and of 
the HI placed according to the halo and particle based model, upper, middle and low row respectively, for the 
CDM (left column) and 1 keV WDM (right column) scenarios. We zoom in a (15 h“^Mpc)^ region, taking a 
slice of 2 /i“^Mpc width. 


we can also infer that most of the hydrogen in the intergalactic medium is ionised, and the neutral 
hydrogen is mainly contained in the so called damped Lyman-a absorber systems (DLAs) with HI 
column densities A^hi above 2 x 10^° atoms cm“^ and secondarily in Lyman limit systems (LLS), at 
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Figure 5. HI column density distribution function, /hi(-^hi), obtained by assigning HI to gas particles em¬ 
ploying two different methods: halo based (dotted lines) and particle based (solid lines) described in 3.1 3.2. 
Observational data are plotted in black for [46] (crosses) and for [47] (circles). We show results at redshifts 
z = 3 (left), z = 4 (middle) and z = 5 (right). In each panel we plot the results for the CDM simulation in 
blue whereas results for the 1 keV WDM cosmology are plotted in red. 


lower HI column densities. In order to test the way we distribute HI in the simulations, we compare 
the abundance of HI absorbers in our simulations with observational data from DLAs [46] and LLS 
[47]. We have assumed that the HI column density distribution function, from observations, does not 
vary in the redshift range studied in this paper. 

The column density distribution function, /hi(-/Vhi )5 is computed for each simulation, redshift, 
and for both HI assignment methods by projecting the HI particles position on a plane and picking 
random lines of sight perpendicular to the plane. For each line the HI column density is computed 
using the SPH radius and HI mass of all gas particles. Details of the computation can be found in 
[31]. 

Results are shown in figure 5 for both the halo based and particle based method, together with 
the observational measurements. For clarity, we only show the results for two cosmological models: 
the models with CDM and 1 keV WDM, since these two are the extreme scenarios and together 
they span the interesting range of haloes that it is also probed by the other models. We find that the 


- 11 - 






























Figure 6. Value of the parameter = Pm{z )/as a function of redshift, obtained by distributing HI 

according to the halo based method (dotted line, all cosmologies) and to the particle based method (colored 
lines). Observational measurements are displayed with error bars in blue, crosses for [48] and triangles for 
[49], 


HI column density distribution obtained by distributing the neutral hydrogen with the two different 
methods considered here is in overall good agreement with the data. 

For the halo based model (dotted lines in figure 5) the agreement works better for column 
densities higher than 10^^ cm“^ while the abundance of absorbers with lower column densities is 
overestimated by this method. We find fhaf our resulfs are redshiff and cosmology dependenf, as fhe 
overesfimafion is weaker af higher redshifls and for low keV mass WDM, which is in fum due fo fhe 
lack of low mass halos in WDM scenarios where low column densify absorbers reside. 

The mismafch in fhe abundance of LLS befween our resulfs and observafions does nof weaken 
fhe resulfs of fhis work, since fheir contribufion fo fhe fofal HI disfribufion is marginal compared fo 
fhe DLAs, which are insfead well reproduced. 

By simulating fhe HI disfribufion using fhe parficle based mefhod (solid lines in figure 5), fhe 
agreemenf is beffer for low HI column densify absorbers while fhe abundance of absorbers wifh large 
column densifies is underestimated af all redshiff and in all cosmologies. 

Generally speaking, fhe fwo differenf Hl-assignmenf mefhods are somehow complemenfary: 
one overesfimafes whaf fhe ofher underesfimafes, leaving fhe dafa poinfs wrapped by fhe fwo ap¬ 
proaches. We sfress fhaf fhe halo based and fhe parficle based mefhods are greatly different both 
in spirit and in the actual details of the implementation as discussed in Sec. 3, but, looking at the 
/hi(A^hi) results, they reasonably well reproduce observations. 

4.2 The HI power spectrum 

We now investigate the impact of WDM on the spatial distribution of neutral hydrogen. We perform 
this task by using the simplest statistical estimator that provides information on the spatial distribution 
of neutral hydrogen: the HI power spectrum. 
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Figure 7. Dimensionless HI power spectrum, Ayj(/c) = fc^PHi(fc)/(27r^), at redshifts z = 3 (left), z = 4 
(middle) and z = 5 (right) obtained by assigning HI to gas particles according to the halo based method. In 
each panel we plot the results of the CDM simulation with a continuous black line and for the WDM ones; 
1 keV mass (dash dotted blue), 2 keV (dashed red), 3 keV (dotted green) and 4 keV (dash dotted magenta). 
The relative difference, in percentage, of the WDM models with respect to the CDM scenario is shown in the 
bottom panels. We set the fc range up to the Nyquist frequency (k ~ 107/iMpc“^). 


We display the HI power spectrum, when the HI distribution is modeled using the halo based 
method, for all the cosmological models studied in this paper in figure 7, while the results when the 
particle based method is used are shown in figure 8. The five different cosmology power spectra 
-Phi(^) get closer to each other, at all redshifts, on large scales, although they never exactly converge 
as happen for the total matter power spectra Pm{k) showed in figure 1. As expected, the lighter 
the WDM mass, the bigger is the discrepancy with the CDM case, with an increase of power. We 
find that whereas the HI power spectrum for the 1 keV model has an amplitude roughly a factor two 
higher than the one from the CDM model, almost independently on the model used to distribute the 
HI, differences between the 4 keV and CDM models are smaller than ~ 10%. 

In the bottom panels of figures 7 and 8, we plot the ratio between the HI power spectra of the 
models with WDM to the one of the model with CDM. In the halo based method, the difference in 
HI power is at least at the ~ 5% level for the 4 keV WDM cosmology at redshift z = 3 and goes up 
to over the ~ 100% for the 1 keV WDM at z = 5. By employing the particle based, the differences 
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Figure 8. Same as figure 7 but modeling the distribution of HI using the particle based method. 


are even higher, reaching more than 300% in the most extreme case of 1 keV at z = 5. 

The reason why the amplitude of the HI power spectrum is higher in the models with WDM 
with respect to the CDM model is straightforward when the neutral hydrogen distribution is modeled 
using the halo based method. By using this method we are forced to distribute the same amount of 
HI in the cosmologies with CDM and WDM. Since we assume that halos with masses smaller than 
Minin do not host HI, and since there is a deficit of low-mass halos in the cosmologies with WDM 
(in comparison to CDM), we need to put more HI into the remaining halos when we assign HI to 
the halos of the WDM cosmologies. Thus, the HI will be more clustered in the models with WDM 
because we have less halos at disposal to be filled in and those available have higher masses and 
therefore larger bias. 

When the HI is modeled using the particle based method it is not that obvious the reason of the 
larger clustering of the neutral hydrogen in the WDM cosmologies. To further investigate this issue 
we have performed the following test. We select all dark matter halos in a given mass interval and 
compute the sum of the HI mass, stellar mass and total mass within those halos for the models with 
CDM and 2 keV WDM. In figure 9 we plot the ratio between the total HI mass and the total mass 
contained in the halos in the selected mass range (top panel), the ratio between the total HI and stellar 
mass (middle panel) and the ratio between the total stellar and total mass (bottom panel). 

Looking at the ratio in the top panel, we see the cutoff of the HI in the high mass 
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Figure 9. Ratios of the total HI mass over the matter Mni/Mtot (upper panel), over the stellar mass 
Mni/Mstars (middle panel) and the total stellar mass over the total mattel Mstars/Aftot (lower panel) cal¬ 
culated as function of the halo mass bins within the ratio is performed, for the CDM (blue lines) and 2 keV 
WDM (red lines) cosmologies with the halo based (dotted lines) and particle based (dashed lines) HI assign¬ 
ment methods at redshift z = 4. 


halos predicted by the halo based model, whereas the particle based results do not show such decline. 
From that panel it is also clear that halos of the same mass have more HI in the models with WDM 
with respect to the CDM model, when the HI is modeled using the halo based method. 
is very similar in the two different models when the HI is modeled using the particle based method. 
These results point out that the amount of HI per total mass in a given dark matter halo is almost the 
same between halos in CDM and WDM cosmologies. Thus, the larger amplitude of the HI power 
spectrum in the models with WDM arises because the bias of the dark matter halos is higher in the 
WDM models. 

The ratio in the middle panel of figure 9 is higher for WDM and comparable 

among the halo and particle based methods for the central mass range bins. The ratio of the stellar 
mass over total mass in the lower panel of course does not change with the different HI 

assignment methods, and it is higher for WDM: it could be sign of more active star formation in the 
WDM halos (but see for example [28]). 

In the next section we analyse in more detail the different bias that HI exhibit in CDM and 


- 15 - 





















10 ° 10 ^ 10 ^ 10 ° 10 ^ 10 ^ 10 ° 10 ^ 10 ^ 

k [h Mpc“^] 


Figure 10. HI bias, 6Hi(fc) = \/ Pm{k)/Pmik), at z = 3 (left column), z = 4 (middle column) and z = 5 
(right column) when the HI distribution is modeled using the halo based (top row) and the particle based 
(middle row) methods. In each panel we plot the results for the COM model with black lines whereas we show 
with colored lines the results for the models with WDM; 1 keV (blue), 2 keV (red), 3 keV (green) and 4 keV 
(magenta). The ratio between the bias of the models with WDM to the model with CDM is displayed in the 
bottom row. We set the k range up to the Nyquist frequency (k ~ 107 h Mpc“^). 


WDM cosmologies, to check whether the total matter clustering properties are indeed reflected by 
the HI clustering properties. 


4.3 The HI bias 

After having computed the total matter and the neutral hydrogen power spectra Pm(^) and Phi(^)> 
respectively, we can now estimate the HI bias, 6 hi(^)> defined as: 




Pm{k) 
Pm{k) ■ 


(4.1) 


In figure 10 we plot the HI bias at redshifts z = 3 (left), z = 4 (middle) and z = 5 (right) when 
the HI distribution is modeled using the halo based (top row) and the particle based (middle row) 
methods. We find that the spatial distribution of neutral hydrogen is more biased in the WDM models, 
with the bias increasing with decreasing WDM mass, at all redshift, for both methods. Among the 
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two different models, the particle based produces a slightly higher HI bias in comparison to the results 
obtained by employing the halo based model. This can be more easily seen in the bottom panels of 
figure 10, where we plot the relative difference in the bias between the models with WDM and CDM, 
for both methods. As discussed in the previous section, the reason why the HI bias is higher in the 
models with WDM is because the neutral hydrogen is more strongly clustered in these models in 
comparison to model with CDM. 


4.4 The 21cm power spectrum and SKAl-LOW forecasts 

Radio telescopes can detect the redshifted 21cm radiation from neutral hydrogen. Therefore, the 
quantity directly measured from observations is not the HI power spectrum, but the 21cm power 
spectrum, which is nothing but the power spectrum of the spatial distribution of neutral hydrogen in 
redshift-space. For concreteness, the 21cm power spectrum is defined as 

P2icm{k) = {snik){6nr{k)), ( 4 . 2 ) 


where 5Tij is given by 

PHl(^ 1 
Phi J ’ 

with phi(^ being the neutral hydrogen density in the redshift-space position s and 


6n{i2) = 6n{z) 


(4.3) 


6Tbiz) = 23.88 XHi 


0.02 


I 0.15 {l + z) 
10 


mK, 


(4.4) 


where xhi = Phi/ph is the average neutral hydrogen fraction. Here we compute the 21cm power 
spectrum for the different cosmological models investigated in this paper and study the sensitivity 
with which the future SKAl-LOW array will be able to discriminate among different models. 

Starting from the spatial distribution of neutral hydrogen in real-space, that we obtain using the 
two methods described in Sec. 3, we then displace the particle positions to obtain the coordinates of 
the particles in redshift-space by doing 


s = af -|- 


l + z 
H{z) 


V\os{x) , 


(4.5) 


with z being the redshift of observation, uios the line of sight component of the peculiar velocity and 
H{z) the Hubble parameter. For convenience, we prefer to work with the dimensionless 21cm power 
spectrum namely: 


^21cm(^) — 


fc^f21cm(fc) 

27r^ 


(4.6) 


We have computed the 21cm power spectrum for each cosmological model and for each redshift 
using each of the two methods employed to model the distribution of HI. For the halo based method 
Dgf = 10-3 is fixed in all scenarios by construction, whereas for the particle based, we have different 
values of for each simulation. 

Since we want to investigate differences between models, and we assume that the value of Dhi is 
fixed by independent observations (such as the abundance of DLAs and LLS), a further normalization 
is needed when computing the 21cm power spectrum, to force all models to have the same value of 
Dhi^: 


(^21cm(^)) 


norm _ a2 /t \ 


Qref 

“HI 

Qsim 

“hi 


(4.7) 


^Notice that by doing this we are assuming that the clustering properties of the neutral hydrogen do not change. 


- 17 - 













where is the value of flni directly obtained by employing the particle based method. 

In figure 11 we show the relative difference in the 21cm power spectrum between the models 
with WDM and the model with CDM. Dotted lines represents the results when the halo based model 
is used whereas solid lines display the results when the particle based model is employed. 

As with the HI power spectrum, we find fhaf fhe amplifude of fhe 21cm power specfrum is higher 
in fhe models wifh WDM fhan in fhe model wifh CDM. Again, fhis is because fhe spafial disfribufion 
of neufral hydrogen is more sfrongly clusfered in fhe models wifh WDM. We have also compufed fhe 
errors wifh which fhe SKAl-LOW radio-felescope will measure fhe 21cm power specfrum. Defails 
of fhe error compufafion can be found in [50]. The shaded areas in figure 11 show fhe quanfify 
where represenfs fhe la error on fhe 21cm dimensionless power 

specfrum of fhe model wifh CDM. Resulfs are shown for observafions limes equal fo to = 1000 hours 
(grey), to = 3000 hours (blue) and to = 5000 hours (fuchsia). 

We find lhal wifh a reasonable observational time of to = 1000 hours [51] the WDM models 
with 1,2, and 3 keV can be distinguished from CDM. On the other hand, the model with 4 keV WDM 
is consistent with CDM at the ~ la confidence level at z = 5. With longer observational times, e.g. 
3000 hours, even the model with 4 keV can be distinguished from the model with CDM at more than 
2a at z = 5. We note that on large scales, the error budget is dominated by cosmic variance, thus, 
the error magnitude barely changes by increasing the observational time. Extending the observation 
time to to = 5000 hours, also the 4 keV case is above the error at ~ 2.5 ct confidence level for 21cm 
emission coming from redshift z = 3. These conclusions are valid for both the methods considered 
in this paper. 

5 Summary and conclusions 

The standard ACDM cosmological model, which has been recently tested with unprecedented accu¬ 
racy by Planck [1], has been proved to be extremely successful in reproducing large scale structure 
observables from the cosmic microwave background to the Lyman-a forest. However, some proper¬ 
ties of the matter distribution at very small scales (< 1 h~^ Mpc) are difficult to reconcile within this 
model. While this putative crisis could be solved by baryonic physics (i.e. feedback), it is intriguing 
to explore other possibilities, like the fact that thermal velocities of the dark matter component could 
affect the properties of matter at small scales (see e.g. the discussion in [9] where other modifications 
of the CDM paradigm are discussed for lifting/solving these problems). The matter power spectrum 
in a cosmology with WDM presents a cut-off on small scales, that arises because WDM can not sig¬ 
nificantly cluster on scales smaller than their mean free path. The tightest constrains on the mass of 
the WDM comes from observations of the Lyman-a forest, which are used to constrain the shape of 
the matter power spectrum on small scales: mwDM ^ 3.3 keV(2cj) [23]. 

In this work we focus on a relatively new observable: intensity mapping. The idea is to perform 
a low angular resolution survey to measure the integrated redshifted 21cm emission from unresolved 
galaxies (see for instance [24]) which will allow to measure the spatial distribution of neutral hydro¬ 
gen, in terms of the 21cm power spectrum. 

In order to do that we investigate the spatial distribution of neutral hydrogen in cosmologies 
with WDM and study the signatures left by those on the shape and amplitude of the 21cm power 
spectrum. We also forecast which with sensitivity will the SKAl-LOW array be able to distinguish 
between different cosmological models. 

Our assumptions rely on the fact that the cosmological parameters are fixed to exquisite preci¬ 
sion by cosmic microwave background experiments, and a further parameter is added that quantifies 
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Figure 11. Relative difference between the 21cm power spectrum of the models with WDM and CDM when 
the HI distribution is modeled using the halo based method (dotted lines) and the particle based method 
(solid lines). Results are shown at z = 3 (left), z = 4 (middle) and z = 5 (right). The error on the 
21cm power spectrum of the model with CDM, normalized to the amplitude of the 21cm power spectrum, 
shown in a shaded region for three different observation times; to = 1000 hours 
(grey), to = 3000 hours (blue) and to = 5000 hours (fuchsia). For clarity, we show the error on 
from one Hl-assignment method only because both are vey similar and overlap at the scale of the plot. 


the suppression of power: we do not explore scenarios in which feedback is important and we as¬ 
sume a reference thermal history that sets the physical properties of the gas in the post-reionization 
era. On one side, our analysis is aggressive in the sense that a more comprehensive analysis of the 
expected signal should probably account for different astrophysical processes (e.g. radiative trans¬ 
fer effects) that could be modeled and marginalized over to quote the final level of significance of a 
WDM model. On the other side, we want to be conservative in two ways: i) we decide to focus on 
a regime which is far from the complex astrophysics of reionization and possibly also not affected 
at an appreciable level by galactic feedback; ii) we use two different models to distribute HI in the 
Universe. Furthermore, this regime has the advantage of being constrained by data (the total amount 
of HI from absorption lines). 

Our simulation suite comprises a set of 5 high-resolution hydrodynamic simulations. We simu¬ 
late a cosmological model with CDM and 4 models with WDM with different particle masses: 1, 2, 
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3 and 4 keV. We model the spatial distribution of neutral hydrogen by assigning, a-posteriori, HI to 
the gas particles in the simulations, by using: i) a halo-based model, in which HI is assigned only to 
gas particles belonging to dark matter halos and ii) a particle-based method in which HI is assigned 
to every single gas particle in the simulation according to its physical properties. We compute the HI 
column density distribution function (CDDF) at redshifts 2 ; = 3,4, 5 and find good agreement with 
observational data in all cases, although models with low WDM masses fail at reproducing the HI 
CDDF for absorbers with low column density when the particle-based model is employed. All mod¬ 
els predict a lower abundance of absorbers with large column densities by using the particle-based 
method. By distributing the HI according to the halo-based model we find an excellenf agreemenf 
wifh observafions for all models, excepf for fhe LLS, whose abundance is over predicfed by using fhis 
mefhod. 

We fhen exfracf fhe HI power specfrum for fhe cosmologies wifh CDM and WDM. Confrary 
fo a naive expecfafion, for which a similar suppression of power should be found in fhe HI and 
21cm power specfra we find insfead thaf fhe amplifude of fhe HI power specfrum is higher in fhe 
models wifh WDM (and increases wifh decreasing WDM mass) fhan in fhe model wifh CDM, af 
all redshifls sfudied in fhis paper and independenfly of fhe mefhod used fo model fhe disfribufion of 
HI. The reason is fhaf in fhe WDM cosmologies we have a deficil of low-mass halos, in comparison 
wifh fhe CDM model (see figure 2), which means fhaf in fhe WDM cosmologies fhe fraclion HI 
residing in fhe mosf massive halos is higher in fhe WDM cosmologies. Therefore, fhe bias of fhe HI, 
^Hi(^) ~ Pm{k)/Pra{k), will increase wifh decreasing WDM mass. This is also in agreemenf wifh 
whaf shown in [52], where fhey infer fhaf WDM low mass halos should exhibif higher baryonic mass 
compared fo fhe CDM case in order fo reproduce fhe local sfellar mass funclion. 

The mosf interesting conclusion is fhaf fhe HI power specfrum of fhe WDM does nof have an 
exponential cul-off, af leasf on fhe scales probes in fhis paper (~ 100 /iMpc“^), fhaf if is however 
presenf in fhe maffer power specfrum. The reason is fhaf neufral hydrogen found resides in halos 
of masses larger fhan ~ 10® and fhe effecf of WDM on fhose halos is fo suppress fheir 

abundance and fhus resulfs in an increased bias. The effecf of WDM on fhe mafter power specfrum 
on small scales is more pronounced since fhe power on fhese scales arises mainly from low mass 
halos. 

Finally, we have compufed fhe 21cm power specfrum for each of our 5 differenf cosmological 
models. As in fhe case of fhe HI power specfrum, we find thaf fhe amplifude of fhe 21cm power 
specfrum is higher in fhe models wifh WDM fhan in fhe models wifh CDM. We find fhaf fhe amplifude 
of fhe 21cm power specfrum is up fo 400% higher for fhe 1 keV case af redshiff z = h and down fo 
10% higher for fhe 4 keV af redshiff z = 3, when comparing if wifh fhe 21cm power specfrum of fhe 
CDM model. 

We have also forecasted fhe observing lime required for fhe SKA 1-LOW array fo be able fo 
distinguish befween models wifh WDM and CDM. Since differences befween fhose models increase 
wifh redshiff, we find fhaf models are more easily disfinguishable af 2 ; = 5 rafher fhan af 2 ; = 3 
independenfly of fhe mefhod used fo model fhe disfribufion of HI. Our resulfs poinl oul fhaf wifh 
1000 hours of observations a WDM model wifh 3 keV can be ruled oul af more fhan 2a at z = 5, 
while SKA 1-LOW will not be able to discriminate among these two models at redshift 2 ; = 3. We 
find fhaf wifh 5000 hours of observafions fhe 4 keV WDM model can be ruled ouf af more fhan lo¬ 
af 2 ; = 3 and af more fhan 2o- af 2 : = 5. In deriving fhese numbers above, we made use only of fhe 
largesl scales k < 1 — 3 h Mpc“^ available, since fhe small scale signal is hindered by noise (see 
figure 11) and Ihereby feedback effecls, if presenf, should affecl fhe signal fo a very small exfenf in 
fhese range of scales and redshifls. 

We slress fhaf in [26] aulhors invesligaled fhe signalures of WDM in fhe 21cm power specfrum 
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during the epoch of reionization. For WDM models with masses of 2 and 4 keV they found differ¬ 
ences in the 21cm power spectrum as large as 30% (for the 4 keV model) and 240% (for the 2 keV 
model), depending on redshift and scale. These differences are much higher that the ones we find in 
this paper since differences between CDM and WDM models increase with redshift. On the other 
hand, in the epoch of reionization the 21cm power spectrum is much severely affected by astrophys- 
ical processes, which should be less important once reionization is over. Thus, observations during 
and after reionization are complementary and will increase the robustness of the constrains on the 
WDM models. 

Regardless of the results obtained for the SKA 1-LOW telescope, our framework explores a 
possible new way to interpret the clustering of matter at small scales and to constrain the fundamental 
nature of dark matter and thereby strengthen the scientific case for conducfing HI/21cm surveys in 
fhe high redshiff universe. 
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